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基于 非 线性 本 构 的 单 壁 矶 纳米 管 的 
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摘 要 :在 近期 石墨 烯 的 力学 实验 中 ,发 现 其 在 有 限 变形 的 条 件 下 应 力 应 变 本 构 关系 是 非 线 性 的 。 
石墨 烯 卷曲 形成 碳 纳米 管 ,大 管 径 单 壁 碳 纳 米 管 的 应 力 应 变 关系 ,在 理论 上 和 石墨 烯 一 致 。 在 本 研 
究 中 ,基于 石墨 烯 的 非 线 性 本 构 关 系 , 建 立 了 新 的 单 壁 碳 纳 米 管 的 Euler-Bernoulli 梁 模 型 。 之 后 对 
两 端 鲜 支 的 情况 ,使 用 Galerkin 方法 和 多 尺度 法 研究 单 壁 碳 纳 米 管 在 均 布 载荷 作用 下 的 静 力 弯 曲 
和 受 连 振动 问题 。 结 果 显 示 , 在 静 力 弯曲 时 ,本 构 中 的 非 线 性 项 对 碳 纳米 管 起 到 了 刚度 软化 的 作 
用 。 在 受 迫 振动 时 ,本 构 中 的 非 线性 项 改变 了 振幅 分 岔 点 的 位 置 。 因 此 本 构 中 的 非 线性 项 对 碳 纳 
=“ 米 管 力学 性 质 的 显著 影响 不 容 忽 视 。 

关键 词 : 单 壁 碳 纳 米 管 ; 非 线 性 本 构 关系 ;Galerkin 方法 ;多 尺度 法 ;Euler-Bernoulli 梁 

中 图 分 类 号 :0322;0369 ”文献 标志 码 :A DOI:10. 11776/j. issn. 1000-4939. 2023.03.019 


Euler-Bernoulli beam model of single-walled carbon nanotubes 
based on the nonlinear constitutive model 


WANG Tengfei ' ,HUANG Kun ,WANG Mingguang' ,GUO Rongxin'” 
(1. Department of Engineering Mechanics, Faculty of Civil Engineering and Mechanics, Kunming 
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Abstract :Recent experiments have shown that the stress-strain relation of graphene is nonlinear under fi- 
nite displacements. Theoretically ,the stress-strain relation of single-walled carbon nanotubes (SWCNTs) 
with big tube diameters are consistent with that of graphene because the SWCNTs are formed by rolling 
graphene. In the present paper, we propose a new Euler-Bernoulli beam model of the SWCNTs based on 
the nonlinear stress-strain relation of graphene. Then ,the static bending and forced vibrations of the SWC- 
NTs are studied by the Galerkin method and the multi-scale method for the case of a hinged-hinged beam. 
The results indicate that the nonlinear terms of stress-strain relation can soften the stiffness of the SWCNTs 
during static bending. In the case of forced vibration, the nonlinear terms of stress-strain relation can 


change the position of the amplitude bifurcation points. Therefore ,the nonlinear term in nonlinear constitu- 
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在 1991 年 ,日 本 名 城 大 学 教授 饭 岛 在 做 实验 
时 ,偶然 发 现 了 一 种 奇特 结构 的 一 维 量子 材料 ,首次 
把 其 命名 为 碳 纳米 管 。 自 此 以 来 ,就 因 其 众多 优异 
的 化 学 和 力学 等 性 能 受到 众多 研究 人 员 的 关注 (9 ， 
应 用 范围 非常 广泛 ,因此 澄清 其 力学 行为 对 于 碳 纳 
米 管 的 应 用 具有 重大 意义 。 碳 纳米 管 力学 行为 相 比 
较 传统 材料 比较 复杂 , 目前 针对 它 的 研究 方法 比较 
有 限 , 主要 分 为 两 种 :理论 分 析 和 实验 研究 。 对 于 实 
验 研究 ,如 TREACY 等 “使 用 TEM 测量 碳 纳米 管 自 
由 端 热 致 振动 的 振幅 , 计算 出 碳 管 的 弹性 模 量 。 
WONG 等 中 在 碳 纳米 管 上 不 同位 置 处 使 用 AFM 探 
钙 诠 加 不 同 作用 力 ,得 到 了 不 同 作用 力 下 的 碳 纳米 
管 财 应 挠 度数 值 。 然 而 ,昂贵 的 实验 设备 和 较 高 的 
实验 失败 率 影响 了 此 种 方法 的 应 用 , 即 在 纳米 尺度 
[SS 前 过 实验 获得 碳 纳米 管 的 力学 性 质 是 极其 困难 。 
轩 必 ,现在 研究 碳 纳米 管 方法 主要 是 理论 分 析 , 主要 
多 法 有 分 子 动力 学 法 、 连 续 介质 力学 等 。 对 于 分 
扰动 力学 方法 ,如 YAKOBSON 等 ”在 研究 碳 纳米 管 
轴 生 弯曲 时 ,采用 数值 计算 方法 -分 子 动力 学 法 得 到 
了 莫 能 量 曲线 及 复杂 的 屈曲 位 型 。KOWAKI 等 
对 不 同 半径 的 碳 纳米 管 在 不 同 温度 下 进行 了 分 子 动 
万 党 模拟 ,并 根据 径 向 分 布 函数 . 均 方位 移 和 原子 构 
型 的 温度 依赖 性 估计 了 “熔化 温度 "。 但 这 种 方法 
计 覆 过 程 复杂 且 耗 时 较 长 ,对 计算 机 运算 能 力 要 求 
极 冬 ,在 成 本 和 时 间 上 受到 很 大 的 制约 ,因此 连续 介 
质 模型 在 纳米 结构 分 析 中 应 用 极其 广泛 。 在 连续 介 
质 力学 中 ,离散 的 原子 结构 被 假定 为 均匀 和 连续 
的 ""。 黄 坤 等 "在 研究 石墨 燃 时 把 其 视 为 均匀 
连续 体 ,并 使 用 薄膜 理论 和 板 过 理论 研究 石墨 烯 的 
力学 性 质 。 黄 坤 等 在 非 局 部 微分 本 构 关系 的 基 
础 上 ,首次 提出 了 具有 小 初始 曲率 的 伯 努 利 - 欧 拉 纳 
米 梁 的 非 线性 偏 微分 积分 方程 模型 。 然 后 应 用 该 模 
型 对 单 壁 碳 纳 米 管 的 静态 弯曲 和 自由 振动 频率 进行 
研究 。 余 阳 等 基于 非 局 部 应 变 梯 度 欧 拉 梁 模型 
研究 了 充 流 单 壁 碳 纳米 管 的 自由 振动 和 波动 特性 。 
ARROYO 等 "中 基于 Cauchy-Born 准则 研究 碳 纳 米 管 
发 现 连 续 介质 模型 与 有 限 元 方法 相 结合 的 模拟 结 
与 涉及 严重 变形 的 零 温度 原子 计算 结果 非常 吻合 。 
YAN 等 ("基于 非 局 部 弹性 理论 推导 了 双 壁 碳 纳 


米 管 的 动力 学 控制 方程 ,并 考虑 了 内 外 壁 之 间 的 范 
德 华 相 互 作用 ,研究 了 压力 驱动 定常 流 的 失 稳 现象 。 

碳 纳 米 管 在 几何 上 可 以 看 成 是 二 维 的 石墨 烯 按 
照 不 同方 向 卷曲 而 成 。 有 研究 表明 , 当 碳 纳米 管 的 
管 径 增 加 达到 一 定 的 数值 时 , 碳 纳米 管 众多 的 力学 
性 能 参数 将 与 石墨 片 层 对 应 的 力学 参数 非常 接近 ， 
比如 弹性 模 量 等 ,因此 石墨 烯 与 碳 纳 米 管 的 应 力 应 
变 关系 理论 上 是 一 致 的 。 在 近期 的 石墨 烯 力学 实验 
中 ,发 现 其 在 有 限 变形 的 条 件 下 应 力 应 变 本 构 关系 
是 非 线 性 的 ”1 。 对 于 单 层 石墨 烯 片 ,LEE 等 …” 通 
过 单 轴 实 验 拉 伸 得 到 其 二 次 非 线 性 本 构 关系 ,之 后 
CADELANO 等 .站 通过 比较 原子 模拟 与 连续 介质 弹 
性 理论 ,得 到 了 LEE 本 构 方 程 中 的 所 有 非 线 性 弹性 
系数 。 本 研究 将 通过 文献 [19 ,21] 中 提出 的 包含 三 
次 非 线性 项 的 非 线性 本 构 关 系 , 来 建立 单 壁 碳 纳米 
管 的 非 线 性 梁 模 型 。 然 后 应 用 该 模型 对 单 壁 碳 纳米 
管 的 静态 弯曲 和 受 迫 振动 进行 了 人 研究 。 


1 单 壁 碳 纳 米 管 的 Euler-Bernoulli 梁 
理论 


本 研究 两 端 不 可 移动 的 单 壁 碳 纳米 管 ( 如 图 1 
所 示 ) 来 建 模 。/ 设 为 矶 纳米 管 的 长 度 ,(x,i) 为 碳 
纳米 管 所 受 的 分 布 载荷 。 


> CS SETS Ces veo 
Wp 

YY x 

SAASA SP SP SP 入 全 了 三? -AP AP 人 


图 1 受 均 布 荷 载 作用 下 的 单 壁 碳 纳 米 管 模型 

Fig.1 Single-walled carbon nanotube model under uniform load 

由 文献 [19 ,21 ] ,可 得 碳 纳米 管 的 非 线性 本 构 为 

0 = 天 ce， + De 十 Ge), (1) 

其 中 :为 二 阶 线 弹 性 系数 ; D 为 三 阶 弹 性 系数 ; G 

为 四 阶 弹性 系数 ; rc 和 ss 分别 为 轴 向 应 力 和 应 变 。 

设 碳 纳米 管 xyz 方 向 的 位 移 分 量 分 别 为 uw、 

w, 这 里 的 v 远 小 于 ww ,因此 本 研究 不 考虑 " 的 影 

响 。 根据 Euler-Bernoulli 梁 理 论 , 可 把 位 移 场 
写 为 ”| 
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加 Ow 和 3 /gwN\ ou 220\ wy) Ou 
Re 0 
(2) 34c (eg) ( 蚂 ) (6) 
其 中 : wu， 、w 和 分别 为 xwywz 方 向 的 位 移 ; :代表 Br) ox 
时 间 。 在 Euler-Bernoull 梁 理 论 中 ,横向 前 切 应 变 和 M=-E, a 加 270[ 至 + sw 下 
横向 法 应 变 将 被 忽略 。 采 用 von Kman 位 移 应 变 pe 
关系 ”1 ,« 方向 的 应 变 为 SA 
_ Nu 9 1 /ow J 90 2 es 
| Ne A 
由 式 (1) 可 以 得 到 有 关 轴 向 力 和 弯 矩 的 表达 Po 本 本 
式 , 分 别 用 字母 NV 和 MM 表示 , 即 A se | 
N = ea4 三 ee, 示 Ds’ + Cea d4 (4) 于 ee Bernoulli 梁 ， 可 以 通过 虚 功 原理 得 到 
M = jea4 = ge ype + Gei)dA (5) a - ma (8) 
A A 多 1 
i; 4 是 碳 纳米 管 的 横 截 面积 。 FM 全 sm 
将 式 (3) 带 入 式 (4) ~ (5) 得 到 轴 力 和 弯 矩 ,为 Ox 和 oF 


La i ee 量 ;(x,t) 是 坐标 轴 的 z 方 
‘SA [+ 二 (并 J+P4[ 兰 + 雪 ( 兰 )] + 向 的 荷载 ,并 忽略 x 方向 上 的 载荷 ,本 研究 不 考虑 
CO 5 we 3 /ou / ow? 方向 上 的 载荷 ,如 图 1 所 示 。 把 式 (6) ~ (7) 代 入 式 

c4[( 闻 ) + 于 (号 ) | a . (8) ~ (9) ,得 " ” 
对 | Ow 9 Ww + D4A|2 dn 2 (2 WY Wt dw 


Ox gx i 2 


OX 


Ox 90x2 OX/ Ox Ox Ox Ox dx’ . 2 ox’ 
3 0 Ou / wd? gu Ow ow 3 w\? 90220 Ou 
G4A|3 (2 aw) a re 王 人 - 2 
[ CA + ( 妆 a 235 ( 汉 tl 2) 3 | lt 
2 gw Fw Ou Ow Ou g220\3 Ow 8220 O31 1 wy)? 82 
> 6LG er A 0 = md! 10 
A ox ox’ Ox Ox ) or’ ? \ox) ox . | of C0) 
证 0w 34w Bu uw su /Pw 1 /ow aw 
+NEE+rF)=- Eh -20D Ey) + + | ) 
Ox Ox’ 和 | 有 dx? Ox Ox Ox gx 9 二 ( 吧 ) Ox 
3 - 346[2 六 a + (2 于 “| - 306|2 (2 | Ow MOLIW, Ou OU Ow 
OX Ox Ox’ Ox’ ox) gx’ Ox ax’ Ox’ Ox 9x2 Ox’ 
Ou Ow 90220 Ow Ow 20 Ow Ow 4 gw Ou Ow Ow 2 
SW 3 ( 中 二 全 3LG Ow 
(区 gx (区 ox’ (0 ox’ pi 上 1 Ox Ox’ a ) 
2 2 2 2 .2 : 3 人 
人 人 -全 人 全 生生 
OX” Ox \ Ox gx OX \ Ax Ox gx \ Ox OX Ox” Ox OX Ox 
1 1 ow] ou 1 /ow 3 ow)? 
EA + DA + C4| |1 王 [ow 
3 /go\ Ou 0220 Ou 3 9220\ 2 / wy? Ow 
2 {9w) 于 a ) Re 7 ( ) } F Sm 11 
4 四 ls g Ox” ox” dx 2 2C Ox” | We or a 
为 了 简化 讨论 , 式 (10) 中 人 铭 去 所 有 的 高 次 非 线 对 于 经 典 细 长 梁 , 纵 向 位 移 v 主要 由 横向 位 移 
性 项 ,得 引起 形变 ,而 且 纵 向 惯性 项 可 以 忽略 。 此 时 ,从 
2 pA es 总 半 2 > 
EA 「2 + +2LD YY = me OE , ， 
O% 0x gx 0Y 07 ot Ou Ow 9 w 9 0 9 2 
Be 2 党 光 2 3 (13) 
(12) Ox OX Ox OX” Ox 
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对 式 (13) 求 一 次 导 , 得 


pe (Se) 
Ox’ x? Ox ox’ 
ow) 0°"w Ow 
( ) =2A LE 4 
人 Ox . Ox Ox 3 
ID 
Ou _ 1 /gw\? 92w)? 1 
| = + c(t) (15) 
» 2 
u =- 记 | ( 妆 ) -A (加 | dx +clCDxX +c (Ci 


(16 ) 
其 中 , ci (CD 和 c(t) 是 时 间 的 困 数 。 对 于 两 疹 贸 文 
的 梁 , 它 的 边界 条 件 为 


u(0,D) = ul,t)=0 (17) 
f pe Pe 
0D = wl,D = a 让 0 =0 
© (18) 
将 过 界 条件 式 (17) ~ (18) 代 入 式 (16) ,可 得 
OO dt) = (19) 


¢ 1 , Ow 2 A | Ow ? 

三 全 ow 全 和 2 
名 ci CD) i | dx + | dx (20) 
将 蒜 (11) 中 高 次 项 全 部 舍 去 ,得 
0340220 人 Ou Ow 


Ou Ow 
i er 3 过 
OX OX Ox OX 


OX Ox 


于 2 | 

3 1 pt 
(总) 上 2 [| 2 
iA 


二 本 有 2 及 4 
304G [2° | + (| | + Fon) = 
Ox \ Ox’ ox) Ox 


(21) 
J SO 


0 w 92w)” Ow 
— +3LG|2 + ) 一 
mm or [2 ls ] 人 ox” Ox 
gw 1 0 F210 0 
12ADD 6ALD {ow) OW 
: ox” > . | Ox 
Ow 0 2W0 Ow 
Bly 一 | -2L,D “|e + F(x,t) 
(22) 
把 式 (20) 带 入 式 (22) 得 


学 学 报 
- 12A7D 5 a -6ALD (到 ) 机 
(人 

= 2 
rN 局 
zp eB 2) s+ A (Es) = Fos, 
(23) 
为 了 方便 下 面 的 讨论 ,对 式 (23 ) 进行 无 量 纲 化 处 
理 , 令 4 =x/1, w=w/1, 7 了 =wot, 其 中 , oo = ee 
m 


这 里 的 w。 表示 两 端 绞 支 梁 的 固有 频率 。 式 (23) 
化 为 
ow 12ALD gw (到 ] 


ox 


6AL,D 人 9 py 和 
OX 


es 276 区 2716 二 
a0 moo Ox mowol dn 


34C | 9°w (3) (3) 0 


mooll ox (ox QO x) ox 


El, gw EA owrl Ow 
一 dx 十 
mol ox | 1 ~” 


mewol 9 x 21 
A nA/ow 1 22D owrl Ow 
中 [ee | + lt 中 a+ 


入 人 人 区) 时] = 局 (24) 
力 界 条件 为 


由 上 可 知 , 该 梁 模 型 边 
_ _ Ow0,7) 9%w(1,7) 
RU 


pm Br 
在 式 (24) 中 是 含有 积分 项 的 非 线 性 偏 微 分 方 
ee 在 此 ,使 用 Calerkin 法 
来 近似 求解 。 在 归 一 化 边界 条 件 下 ,可 设 微分 方程 
式 (24) 的 近似 解 ” 为 


Dy nSinnm% (25) 


里 只 取 第 一 项 ,将 式 (25) 带 入 式 (24) 中 ， 再 在 方 
程 的 两 边 同 时 乘 snm xz, 然后 在 区 间 [0,1] 上 积 4 
( Galerkin 一 次 截断 ) ,可 得 到 描述 单 壁 碳 纳 米 管 一 
阶 模 态 受 迫 振动 的 微分 方程 ,为 


gn, Eom EAm’ | a 31Cm 
9 mool 4mo2l? 忆 4meaw2ls 
万 DT 
1 3 = rR (x,t)sinmxdx (26) 
2mowol mool 
上 式 可 写 为 

D+hn+kn = 下 (27) 
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其 中 
EAm’ 4Am 3LGmw ALD™’ 
3 一 sa(1 + 2 bb 276 一 276) 
mwo l 4mwo ! 2mowo 1/ 
Elm™’ 
ke? = 0 2 | (x,1)sinmxdx 
mowo l mowo l 


2 静 力 学 分 析 


去 掉 式 (27 ) 的 惯性 项 ,可 得 碳 纳 米 管 的 静 力 学 
方程 
En+knm =F (28) 
在 此 为 方便 讨论 ,本 研究 以 (15 ,15 ) 碳 纳米 管 
为 例 ,d =2.034 nm 为 其 直径 的 大 小 ,i, =0.34 nm 为 
人 它 的 物理 参数 为 : E = 1 TPa ; G = 
6. 58 TPa [2 。 其 它 参数 为 :1 = 11 nm;m = 4.866 x 
16 Wkg/m;l, = 1.154 nm ;7 = 0.953 nm’;A = 
20191 nm ;wo = 1.5778 x10”。 
4 和 2 些 参数 代入 得 到 
k, =1,1 =56.8826,D =G=0 
k, = 1,k = 12.3536,D#0,G=0 
k, = 1,k, = 15.6344,D#0,G0 


通过 式 (28) 和 上 述 的 数值 , 即 可 得 到 碳 纳米 管 
1 人力 寺 遇 扳 朵 如 下 图 2 所 示 。 
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图 2 碳 纳米 管 静 力 弯曲 变形 
Fig.2 Bending deformations under static loads 
由 图 2 可 知 , 碳 纳 米 管 应 力 应 变 非 线 性 本 构 中 
的 二 次 和 三 次 非 线性 项 均 对 碳 纳 米 管 的 静 力 学 性 质 
有 具有 非常 明显 的 影响 ,这 种 影响 程度 会 随 外 激励 的 
增加 而 加 深 。 


3 弯曲 振动 分 析 


在 式 (26) 加 阻尼 项 ,可 得 碳 纳米 管 在 一 


阶 模 态 
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下 的 有 阻尼 受 迫 振动 方程 , 即 


El, mr 4 2 
on Ey 2 + | 2 (+ 2)+ 
9 ot mowo’l 4moo 大 l 
3LGm AD |] ;3 


Se 2 | (xX,t)sinmxdx 
4mowo LD 2mowo 


mowol 
(29) 
其 中 , C 为 阻尼 常数 。 
为 了 便于 摄 动 求 解 方程 , 令 C = 2ec， 
(Xx,1) = F(x)coswt 。 因 此 上 式 方程 可 写 为 


e+ = efeoswi (30) 


oF 
其 中 
EAm™’ 4AT2 3lhGm™ ALD™’ 
= 十 5 ; 
Im 忆 ) 4mowils 2mw2l 
EI 
hi = 2 人 .六 (x)sinmxdx; 
mool 


Pe wg 为 阻尼 和 常数。 在 此 
使 用 多 尺度 方法 求解 式 (30) ,引入 频 差 r , 令 ec = 
w 一 ,把 方程 的 解 在 两 个 时 间 尺 度 上 展开 为 
n I (31) 
其 中 , 7, = ei,(n = 0,1.2，… 
se es ,2 的 系数 相等 
得 到 和 | 


1 


ea : Don +hn =0 (32a) 
es: Do + hn(t) = -2D,D,n, -2cDono - 

kn +feos (kT + oT,) (32b) 
其 中 , Do,D; 分 别 表示 对 7 ,7 求 导 。 式 (32a) 的 解 
可 表示 为 


mo = A(T,)exp(ikiTo) + ee (33) 
其 中 , cc 表示 其 左边 各 项 的 共 红 复数 。 
把 mo 代入 式 (32b) ,将 cos(h T+o7T,) 以 指 
数 函 数 形式 表示 ,得 到 
- [2ik (DA + cA) +31424] ， 
exp(ikT) -kA iexp(3ikT ) + 


Dem, + kin S 


Ffexp[i(hiT, ps 


其 中 , cc 表示 其 左边 各 项 的 共 恩 复数 。 
为 消除 式 (34) 中 的 永年 项 ,要 求 函 数 4 满足 
2ik (CD,A + cA)+ 3hsA’A - 3fexp liaT,) =0 


(35) 
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将 复 函 数 4 表示 为 指数 形式 
4 = 二 aexp(ig) (36 ) 


2 
其 中 ,a 和 8B 是 待定 的 关于 7, 的 实 函 数 。 把 所 得 结 
果 的 实 部 和 虚 部 分 开 , 得 到 a 和 一 阶 常 微分 方 


程 组 


1f 
Da = -ca + a sin(g7T, - B) (37a) 
aD,B = 3 as = 上 7T, -B) (37b) 
2 8 hk 2 hk OIL Bb 


将 式 (36) 代 入 式 (33), 则 解 式 (31) 得 第 一 
近似 解 


7 = eacos(kt +B)+ O(a’) (38) 
其 中 ,0( ) 表 示 为 同 阶 无 穷 小 。 


| 工 了 ， 
Da = 一 ca + 7 hy (39a) 
aDy = oa + 了 = ep (39b) 

8 hk 2 Ek 
令 Dya = D,y = 0 ,系统 的 稳 态 解 为 
_ 工 了 ， 

co = 了 (40a) 

3 局 ， 二 广 
8 EE ag = 一 和 (40b ) 


上 上 击 二 式 的 平方 后 求 和 得 
Se 


站 4 
设立 # 0 时 ,把 式 (41) 整 理 为 
三 士 一 42 
+ (ja °) 2 


这 里 仍 以 (15 ,15 ) 单 壁 碳 纳米 管 为 例 , 基 本 参 
数 与 静 力 分 析 所 举 实例 一 致 。 对 于 碳 纳米 管 的 阻尼 
系数 ,为 了 简化 讨论 ,在 此 取 C = 0.05 ,es =0.1， 
f=5 得 到 对 应 的 频 响 曲 线 如 图 3 所 示 。 

由 图 3 可 知 , 当 o 处 于 某 一 临界 值 时 ,振幅 将 出 
现 分 岔 ,而 本 构 中 的 非 线性 项 改变 了 振幅 分 贫 点 的 
位 置 。 在 分 岔 点 处 , 碳 纳 米 管 的 振幅 会 伴随 参数 的 
轻微 变化 而 突然 改变 ,这 对 结构 振动 产生 了 显著 影 
响 。 为 考查 摄 动 解 的 准确 性 ,对 式 (29 ) 进行 数值 计 
算 。 在 图 3 上 选取 点 o = 3.76 ,对 3 种 不 同 参数 进 
行 数值 计算 。 结 果 显 示 , 摄 动 解 准确 地 描述 碳 管 的 
受 迫 振动 行为 ,如 图 4 所 示 。 此 外 ,图 3 还 显示 ,本 
构 中 的 非 线性 项 起 到 刚度 软化 的 作用 。 


(3.76，0.98) 
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图 3 外 激励 1 = 5 时 的 频 响 曲线 
Fig.3 Frequency response curve forf = 5 
—— D=0, G=0 
0.12| —a— D0, G-0 
—@— D@#0, GA0 
0.09 i 
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图 4 相位 图 在 o = 3.76 
Fig.4 Phase diagram foro = 3.76 
通过 式 (41) ,并 取 定 o = 5 的 数值 ,可 得 到 对 应 
的 激励 振幅 响应 曲线 ,如 图 5 所 示 。 
1.4 (7.32，1.200) 


一 D=0, G=0 
12| —a— DA#0,G=0 


—e— DzA0, G0 


图 5 当 c = 5 时 的 外 激励 响应 曲线 

Fig.5 Response curve of external excitation foror = 5 

该 图 显示 ,对 于 相同 的 er 值 , 随 着 外 激励 值 的 增 
加 , 啊 应 曲线 都 出 现 了 分 岔 点 ,但 分 岔 点 的 位 置 并 不 
同 。 这 说 明 本 构 中 的 非 线性 项 同样 改变 了 这 种 情况 
下 分 岔 点 的 位 置 。 为 检验 解析 解 的 准确 性 ,在 图 5 
上 选取 f = 7. 32 对 不 同 的 参数 对 式 (29 ) 进行 数值 计 
算 ,结果 同样 显示 ,解析 解 具 有 良好 的 精度 ,如 网 6 
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图 6 当 o = 5 时 的 相位 图 
Fig.6 Phase diagram foro = 5 


论 


Dr 


中 在 本 研究 中 ,我 们 通过 石墨 烯 的 非 线 性 应 力 应 
或 奖 系 和 经 典 Euler-Bernoulli 梁 的 平衡 方程 ,建立 了 
新 的 单 壁 碳 纳米 管 的 非 线性 梁 模 型 。 之 后 通过 
rkin 方法 和 多 尺度 法 研究 了 两 端 饺 支 条件 下 ， 
(315) 单 壁 碳 纳 米 管 在 垂直 于 轴线 的 分 布 力作 用 
的 静 力 弯曲 和 一 阶 模 态 受 迫 振动 问题 。 本 研究 得 
到 NI 点 主要 结论 。 
"1 ) 本 构 中 的 非 线 性 项 对 碳 纳米 管 的 静 力 学 有 显 
车 的 影响 ,起 到 了 刚度 软化 的 作用 , 随 着 外 载荷 的 增 
-2 ) 在 受 迫 振 动 时 ,本 构 中 的 非 线性 项 改变 了 振 
贫 点 的 位 置 , 这 对 结构 振动 将 产生 显著 影响 。 
“3) 本 研究 结果 表明 ,对 基于 线性 本 构 关系 的 单 
壁 碳 纳 米 管 研究 ,有 必要 重新 检视 其 准确 性 。 
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